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Abstract 

Exact solutions are derived for an n-dimensional radial wave equation with a gen- 
eral power nonlinearity. The method, which is applicable more generally to other 
nonlinear PDEs, involves an ansatz technique to solve a first-order PDE system of 
group-invariant variables given by group foliations of the wave equation, using the 
one-dimensional admitted point symmetry groups. (These groups comprise scalings 
and time translations, admitted for any nonlinearity power, in addition to space- 
time inversions admitted for a particular conformal nonlinearity power). This is 
shown to yield not only group-invariant solutions as derived by standard symmetry 
reduction, but also other exact solutions of a more general form. In particular, solu- 
tions with interesting analytical behavior connected with blow ups as well as static 
monopoles are obtained. 
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1 Introduction 



Exact solutions are of considerable interest in the analysis of nonlinear par- 
tial differential equations (PDEs) - - in particular, similarity solutions and 
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more general group-invariant solutions typically provide insight into asymp- 
totic behavior, critical or blow-up behavior, and also give a means by which 
to test numerical solution methods. For a nonlinear PDE in two independent 
variables, invariance under a one-parameter point symmetry group leads to 
a reduction of order to a nonlinear ordinary differential equation (ODE) for 
the corresponding group-invariant solutions [1,2]. The strong utility of this 
method lies in the fact that the point symmetries admitted by a given nonlin- 
ear PDE can be found, often in explicit form, by an algorithmic procedure (and 
nowadays there are many symbolic software programs available for computa- 
tion of point symmetries). As a consequence, all point-symmetry reductions 
to nonlinear ODEs are readily derived for a given nonlinear PDE. However, it 
remains necessary to solve these ODEs to obtain the group-invariant solutions 
explicitly [3] and in most cases this is still very difficult. Indeed, in general, it 
may be necessary to resort to special ansatzes or ad hoc solution techniques 
in order to find explicit solutions. 

In this paper we explore an alternative method that is founded on group- 
invariant foliation equations [4] to find exact solutions of nonlinear PDEs. To 
set up the method, for any admitted point symmetry, the first step involves 
converting the given nonlinear PDE into an equivalent first-order PDE sys- 
tem whose independent and dependent variables are given by, respectively, the 
classical invariants and differential invariants of the point symmetry transfor- 
mation. This yields what is known as a group-resolving system, since formally 
it amounts to constructing a group foliation of the nonlinear PDE by employ- 
ing jet space coordinates adapted to the given point symmetry group. Viewed 
geometrically, this system quotients out the action of the given point symme- 
try group, so the points in its solution jet space correspond (generally in a 
one-to-one fashion) to the orbits of the symmetry group in the solution jet 
space of the given nonlinear PDE [5]. Each orbit is determined by integra- 
tion of the system of defining equations for the differential invariants of this 
symmetry group. These equations have a well known automorphic property 
[6] that all of their solutions are related through a transitive group action by 
symmetry transformations, allowing the equations to be explicitly integrated 
once any particular solution is known. Hence the problem of finding exact 
solutions of the original nonlinear PDE is reduced just to seeking particular 
explicit solutions of the group-resolving system. As the next step, the key 
idea of our method is to solve the group-resolving system by a separation of 
variables ansatz technique that is tailored to the form of nonlinearities in the 
PDEs of this system, expressed in terms of the group-invariant variables. In 
practice, with this technique, many explicit solutions are easily found, whose 
form would not be readily obvious just by trying simple ansatzes using the 
original independent and dependent variables in the nonlinear PDE, or by 
simply writing down the form for classical group-invariant solutions. Most 
important, because all solutions of the nonlinear PDE are contained in each 
of its group-resolving systems, the method is not limited to obtaining exact 
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solutions that just arise from classical point-symmetry reduction. 

We will apply this systematic method to a physically and analytically inter- 
esting n-dimensional radial wave equation of a semilinear form with a general 
power nonlinearity. In Sec. 2, the point symmetries and corresponding group- 
invariant reductions for this wave equation are considered. We find that the 
nonlinear ODEs given by the reductions admit extra point symmetries in 
some cases, leading to reduction to quadrature for certain group-invariant so- 
lutions of the wave equation. These solutions in particular include finite energy 
monopoles and static "solitons", which are of analytical interest. In Sec. 3, 
the method using group-resolving systems associated with the admitted point 
symmetries of the wave equation is introduced, employing our ansatz tech- 
nique. Exact solutions obtained from this method are summarized in Sec. 4, 
along with a discussion of their analytical features of interest. Our results give 
several new families of time-dependent solutions for n > 1, including some 
with analytical properties that pertain to blow-up phenomena in n > 2 di- 
mensions. (Additional results for the case of n = 1 dimensions will be reported 
elsewhere.) Finally, some concluding remarks are made in Sec. 5. 



2 Point-symmetry reductions 

We consider the semilinear radial wave equation —u tt + u rr + ' l] -^-u r = ±u\u\ q ^ 1 
for u(t,r), with a nonlinear interaction described by a power q ^ 1, where r 
denotes the radial coordinate in n > 1 dimensions. Global (long-time) behavior 
of solutions to the Cauchy problem for this wave equation has attracted much 
analytical interest over the past few decades (see Refs. [7,8] for a summary). 
Note that stability of solutions relies on the conserved energy S = j™{\u t 2 + 
|w r 2 =t ^-|M| 9+1 )^ n_1 dr to be positive definite, which requires the sign of the 
nonlinear interaction term in the wave equation to be positive. In this case, 
it is known that any smooth initial data u(0, r), Ut(0, r) evolves to a smooth 
global solution u(t,r) for all t > provided the nonlinearity power q is less 
than a critical value 



(in particular this includes all powers q < oo when n — 2). The same result has 
been proved to hold in the critical case g* = 5 when n = 3. For supercritical 
powers q > q* when n > 3, it is widely expected based on numerical studies 
that some smooth initial data will evolve to a solution u(t,r) that blows up 
at a finite time t = T < oo. In light of these results, exact solutions of this 
wave equation may be of definite analytical and numerical interest. 



q* = (n + 2)/(n - 2) = 1 + 4/(n - 2) 
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For our purpose it will be useful to work with a slightly modified nonlinear 
interaction term, 

u tt — u rr — (n — l)r _1 w r = ku q , k = ±1, (2.2) 



which is easier to handle for finding exact solutions. Note, here, we will regard u 
to be a real- valued function, as is compatible with the case of all integer values 
q; in the case of non- integer values q, we will either restrict the domain of u(t, r) 
so that u > 0, or allow a formal complexification of u(t, r). In any case, it will 
be straightforward to patch together such solutions, if necessary, to satisfy the 
more analytically well-behaved wave equation —u tt + u rr + IL ^-w f , = ±u\u\ g ~ 1 . 

The geometrical point symmetries admitted by the wave equation (2.2) for 
n > 1 are well known [7]: 



time translation X trans . —d/dt for all q, (2.3) 
scaling X scal . = td/dt + rd/dr + -^—ud/du for all q ^ 1, (2.4) 

space-time inversion Xi nver . = (t 2 + r 2 )d/dt + 2rtd/dr + (1 — n)tud/du 

only for q = (n + 3)/(n — 1), (2.5) 

where X. is the infinitesimal generator of a one-parameter group of point 
transformations acting on (t,r,u). The special power for which the inversion 
exists is commonly called the conformal power 

q c = (n + 3)/(n- 1) = 1 + 4/(n - 1). (2.6) 



To begin, we discuss the group-invariant solutions arising from reduction of 
the wave equation (2.2) to a nonlinear ODE under each of these point symme- 
tries. It will be worth pointing out that these ODEs each inherit a variational 
structure related to the natural action principle 

+oo OO 

S[u] =11 {-u t 2 + u r 2 — u q+l )r n ~ l dr dt (2.7) 

J J Q ~\~ X 

-oo 



for this wave equation. Note that, as we easily see, this action principle is 
formally invariant under both the time translation and inversion symmetries, 
for any dimension n, but not under the scaling symmetry unless the dimension 
is n — 1 + 4/(q — 1) corresponding to the conformal power q = q c . 
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2.1 Time translation invariance 



The form for group-invariant solutions is obviously u = U(r), satisfying the 
ODE 

U" + {n- l)r- l U' + kU q = 0. (2.8) 

Despite its simplicity, this nonlinear second-order ODE is intractable to solve 
explicitly in general. However, for special values of the parameters n and q, 
it will be possible to integrate (2.8) completely to quadrature if it admits a 
variational point symmetry (see Ref. [1] for details), i.e. a symmetry of ODE 
(2.8) that leaves invariant the reduced action principle 

/9k 
(U' 2 - dr (2.9) 

^ 

for this ODE. The computation of point symmetries is straightforward and 
leads to the following result. 

Proposition 1: For translation-invariant solutions of wave equation (2.2), the 
point symmetries admitted by ODE (2.8) for n > 1 are comprised by 

scaling X scal = rd/dr + — Ud/dU for all q ^ 1, (2.10) 

1-q 

and 

non-rigid dilation X dn . = r 2 ~ n (rd/dr + (n - 2)Ud/dU) 

only for g = (4 - ra)/(n - 2), n ^ 2. (2.11) 

The dilation is a variational symmetry for all n 7^ 2 (i.e. such that q 7^ 00), 
while the scaling is only a variational symmetry when n = 2(q+l)/(q — 1) 7^ 2 
corresponding to the critical power q = 1 + 4/(n — 2) = q*. 

The appearance of special powers, especially the critical power, is interesting 
because these powers are not distinguished by the geometrical point symme- 
tries of the wave equation. Note that the non-rigid dilation generalizes a uni- 
form dilation group, r — > \r, u — > flu, such that A, Q are no longer constants 
but depend on a power of r. 

To carry out the reduction to quadrature implied by Proposition 1, we change 
variables to canonical coordinates derived from the coefficients of the in- 
finitesimal transformations X. (see Ref. [2] for details) in the variational case. 
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For the scaling symmetry (2.10), these coordinates are given by x = lnr, 
v = r~ 1+n / 2 U, yielding the following ODE for v(x) 

v" - p 2 v + kv q = 0, p= l-ra/2, q = l + 4/(ra-2) (2.12) 

(i.e. a nonlinear oscillator equation) which is easily integrated to quadrature. 
Thus, we have the implicit general solution (n ^ 2) 

x + c = ± J l/yj2c + (1 - n/2) 2 v 2 - k(l - 2/n)v 2n /( n ~ 2 )dv (2.13) 

with arbitrary constants c, c. This includes as a particular case the invariant 
solution v = const of reduced ODE (2.12) under the scaling symmetry (2.10), 
given explicitly by 

v = ((l-n/2) /Vk) n/2 -\ (2.14) 

Canonical coordinates for the dilation symmetry (2.11) are given by x = 
_L_ r n-2^ y _ r n-2jj _ yields a similar integrable ODE for v(x) 

v » + k v i = 0, g = -l + 2/(n-2) (2.15) 

which has the implicit general solution (n ^ 2) 

x + c = ± y l/y^c - k(n - 2)v 2 l^- 2 Uv (2.16) 

with arbitrary constants c, c. We note, here, the invariant solution of reduced 
ODE (2.15) under the dilation symmetry (2.11) is trivial, v — 0. 

Thus, symmetry reduction yields two translation-invariant solutions, in im- 
plicit general form (2.13) and (2.16), for special powers q — l + 4/(n — 2) = g* 
and q = — 1 + 2/(n — 2), respectively. In the case of other powers, although 
the ODE (2.8) admits the scaling symmetry (2.10) for all q ^ 1, and therefore 
its order can be reduced through a symmetry reduction procedure [1,2], the 
reduced first-order ODE is not tractable to solve by any standard integrating 
factor methods [2], apart from the invariant case analogous to the solution 
(2.14). Consequently, it would appear to be necessary to try special ansatzes 
or ad hoc techniques to solve the second-order ODE (2.8) to obtain additional 
translation-invariant solutions. 

A similar situation occurs for the other group-invariant solutions, as we will 
now see. 
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2.2 Scaling invariance 



The form for group- invariant solutions is easily seen to be u = t p £/(£), where 
p — 2/(1 — q) and £ = r/t, with [/(£) satisfying the nonlinear second-order 
ODE 

(1 - f)U" + ((n - l)r' + 2(p - 1)0^' + p(l - p)^ + = 0. (2.17) 



This ODE is considerably more complicated than the one for translation- 
invariant solutions. Nevertheless, it will still be possible to reduce (2.17) to 
quadrature for special values of the parameters n and q if it admits a varia- 
tional point symmetry. An added complication here is that the action prin- 
ciple (2.7) does not have a well-defined scaling reduction (due to its lack of 
invariance under the scaling symmetry), but this difficulty can be by-passed 
if we consider multipliers h(£) for the ODE (2.17) (see Ref. [2]). We find that 
h(£) = (1 — £ 2 )-p+( 1 -™)/ 2 leads to the existence of an action principle 

+ r° 2k 

s sca ,[u}= ( P (i - P )u 2 - (i - ew 2 + — -t/^xor- 1 ^- (2.i8) 

J q + 1 

— oo 



Now, a straightforward computation of point symmetries gives the following 
result. 

Proposition 2: For scaling-invariant solutions of wave equation (2.2), the point 
symmetries admitted by ODE (2.17) consist of only 



non-rigid dilation X dil . = J\l - £ 2 |(£d/d£ + (1 - n/2)Ud/dU) 

for q = (ra + 2)/(n-2), n ^ 2. (2.19) 

This symmetry is variational for n = 2(q + l)/(q — 1) ^ 2 corresponding to 
the critical power q = 1 + 4/(n — 2) = for all n^2. 

Thus, ODE (2.17) has a reduction to quadrature precisely in the case of the 
critical power. Changing variables to canonical coordinates, given by v — 
t-i+n/2-fj = iin^v^ if iti < i an d x = e _ arctan i = if |£| > f 

(with x smooth at £ = 1), we obtain the following ODE for v(x) 

sv" - p 2 v + kv q = 0, p=l-n/2,g=l + 4/(n-2) (2.20) 
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where s = sgn(l — £ 2 ). This nonlinear ODE has the same form as ODE (2.12), 
with the implicit general solution (n 7^ 2) 

x + c = ± J l/y^c + s(l - ra/2)V - sk(l - 2/n)v 2 "/( n - 2 )cfc (2.21) 

with arbitrary constants c, c. As a special case, the invariant solution v = const 
of reduced ODE (2.20) under the dilation symmetry (2.19) is given by the same 
expression (2.14) as for ODE (2.12). 

Hence, other than the special case (2.14) and its obvious generalization to non- 
critical powers, symmetry reduction yields just the scaling-invariant solution 
in implicit general form (2.21), for the critical power q = 1 + 4/(n — 2) = q*. 

2.3 Conformal invariance 

The form for group-invariant solutions here is given by u = r p £/(£) with p = 
2/(1 — q) = (1 — n)/2 and q = q c = 1 + 4/(n — 1) where £ = (t 2 — r 2 )/r is the 
invariant associated with inversions (radial conformal transformations) on the 
space and time variables. The ODE satisfied by U(£) is readily found to be 

£ 2 U" + 2£U' l)U + kU q = $. (2.22) 

In contrast with the scaling case, it is straightforward to derive an action 
principle 

+ T 2k 
S illvc ,[U] = / {p{p + l)U 2 + eU' 2 - -^jU" +1 ) d£ (2.23) 

J q + 1 

—00 

for the nonlinear second-order ODE (2.22), without need for a multiplier, due 
to the dilation invariance of the action principle (2.7). Proceeding as before, 
we now determine the variational point symmetries admitted by this ODE. 

Proposition 3: For conformally-invariant solutions of wave equation (2.2), the 
only point symmetry admitted by ODE (2.22) for n > 1 is a non-variational 
scaling of its independent variable 

scaling X scal . = id/di. (2.24) 

Because this scaling is not a variational symmetry, we do not obtain a complete 
reduction of the ODE (2.22) to quadrature. In particular, with a change of 
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variables to canonical coordinates x — m£, v — U, the ODE becomes 

v" + v' -p{p + l)v + kv q = 0, p=(l-n)/2,g = l + 4/(n-l) (2.25) 

(i.e. a damped nonlinear oscillator equation). This ODE is intractable to solve 
explicitly in general, except in the case of the invariant solution v = const 
under the scaling symmetry (2.24), which yields 



However, without resorting to special ansatzes or other techniques, we do not 
obtain any other explicit conformally-invariant solutions through symmetry 
reduction. 

2.4 Group action on solutions and invariance under optimal subgroups 

We now discuss some deeper group-theoretical aspects of our results on in- 
variant solutions of wave equation (2.2). First it is natural to look at how 
these solutions transform under the full symmetry group, Q q , generated by 
the geometrical point symmetries (2.3), (2.4), (2.5). Note the group struc- 
ture depends on the nonlinearity power q ^ 1. The one-dimensional point 
symmetry subgroups of translations, scalings, and inversions are shown in the 
Appendix. 

For all q ^ 1 other than the conformal power q c , the point symmetry algebra 
is given by the commutator of the translation and scaling generators 

[X trans ., X sca J = X trans .. (2.27) 

This generates a solvable two-dimensional symmetry group Q q ^ qr — U(l) x 
U(l) consisting of a semidirect product of the scaling subgroup (A. 2) and the 
translation subgroup (A.l). The explicit families of translation-invariant so- 
lutions (2.13) and (2.16) found for special powers of q each remain invariant 
under the action of this group, subject to a shift (c, c) — > (c + In A, c) and 
a scaling (c, c) — > (c\ 2 ~ n , cA^ 4-2 ™^*™ -3 )) on the integration constants. In con- 
trast, the explicit family of scaling-invariant solutions (2.21) is non- invariant, 
as it gives rise to a one-dimensional orbit under the time translation subgroup, 
with the integration constants remaining unchanged. 

For the case of the conformal power q = q c , the point symmetry algebra is 
enlarged by commutators involving the inversion generator, 

[Xt rans ., X; nvcr J 2X sca i , [X sca j Xi nvcr J Xi nvcri . (2.28) 



v = 



((n-l)(n-3)/(2v / A?)) (n ~ 1)/4 . 
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The commutators (2.27) and (2.28) together generate a semisimple three- 
dimensional symmetry group G q = qc — SX(2,]R) obtained by composition of 
the inversion subgroup (A. 3) with the scaling and translation subgroups (A. 2) 
and (A.l). Under the action of the full group, the single explicit solution (2.26) 
found for the conformal power q c is invariant. 

An interesting feature arises from the three-dimensional nature of the symme- 
try group structure in the conformal power case. With respect to conjugation 
by the full group, the translation and inversion subgroups fall into the same 
conjugacy class, while the scaling subgroup comprises its own conjugacy class. 
The origin of this structure is easily explained by means of the following special 
conformal point symmetry of wave equation (2.2), 

u(t,r)^(t 2 -r 2 )- p u(-t/(t 2 -r 2 ),r/(t 2 -r 2 )), p = 2/(1 - q c ) (2.29) 

which corresponds to an inner automorphism of the symmetry group SL(2, R), 
as shown in the Appendix. This transformation acts as an involution X trans <-> 
— X invcr , X sca i <-> — X sca i , on the symmetry algebra. Thus it follows that a 
complete enumeration of conjugacy classes of the symmetry group is rep- 
resented by the one-dimensional subalgebras X tra ns. or X inver ., X sca i., and 
Xtrans. + X inver . These are referred to as the optimal subalgebras for group- 
invariance considerations [9,1]. 

As a consequence, through the point symmetry transformation (2.29), ODE 
(2.8) for translation-invariant solutions with q = q c is equivalent to ODE (2.22) 
for conformally-invariant solutions. An essentially different ODE will arise for 
solutions invariant under the subgroup generated by 

X trans . + X invcr . = (1 + t 2 + r 2 )d/dt + 2rtd/dr + (1 - n)tud/du (2.30) 

combining translations and inversions, as we now investigate. 

The form for such group-invariant solutions is a generalization of conformally- 
invariant solutions, given by u = r p £/(£) where, now, £ = (1 + t 2 — r 2 )/r, with 
p — 2/(1 — q) = 1 — n/2 and q — q c — 1 + 4/(n — 1). The corresponding 
reduction of wave equation (2.2) yields the ODE 

(£ 2 + A)U" + 2£U' - p{p + l)U + kU q = 0. (2.31) 

Compared to the conformal case, scaling invariance is obviously lost in this 
reduction. A computation of point symmetries for this ODE shows that the 
extra dilation invariance holding in the translation case is also lost. 

Proposition 4: ODE (2.31) for group-invariant solutions of wave equation (2.2) 
admits no point symmetries for n > 1. 
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This result is perhaps unsurprising because of the less geometric nature of 
the symmetry group generator (2.30) here. The lack of invariance of ODE 
(2.31) means it has no symmetry reductions. It is also intractable to solve by 
integrating factor methods [2], and hence we do not find any explicit solutions. 



3 Group foliation method and ansatzes 



The basis of our method is the introduction of group foliation equations associ- 
ated with the point symmetries admitted by the wave equation (2.2). There is 
an algorithmic construction for group foliations [4] in terms of group-invariant 
jet space coordinates given by the invariants and differential invariants of an 
admitted one-dimensional point symmetry group. To explain the ideas, we 
first consider the scaling symmetry (2.4). 

A convenient choice of invariants and differential invariants is given by 

x = t/r, v = r~ p u (3.1) 



and 



G 



H 



r l ~ p u r 



(3.2) 



where p = 2/(1 — q). While x and v here are mutually independent, G and H 
are related by equality of mixed r,t derivatives on u t and u r , which gives 

D r {r p ^G) = D t (r p - l H) (3.3) 



where D r , D t denote total derivatives with respect to r, t. Furthermore, v, G, H 
are necessarily related through the wave equation (2.2) by 

Dtir^G) - D r (r p - l H) = r p ~ 2 ((n - l)H + kv q ). (3.4) 



Now we put G = G(x, v ), H = H(x, v) into equations (3.3) and (3.4) to arrive 
at a first-order PDE system 



(p - l)G - xG x + (H- pv)G v -H x - GH V = 0, (3.5) 
G x + GG V - (p + n - 2)H + xH x - (H - pv)H v = kv q , (3.6) 

with independent variables x,v, and dependent variables G, H. These PDEs 
will be called the scaling-group resolving system for the wave equation (2.2). 
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We establish, first of all, that this system (3.5)-(3.6) contains all scaling- 
invariant solutions of the wave equation (2.2). Consider any solution 



u(t,r) = r p U( y x), (3.7) 

with scaling-invariant form, where 

(1 - x 2 )U" + (2p + n- 3)xU' -p(p + n- 2)U - kU q = (3.8) 

is the ODE given by reduction of PDE (2.2). From relation (3.2) we have 

G — U', H = P U - xU', (3.9) 

which together yield 

H = pv-xG. (3.10) 



This relation is easily verified to satisfy PDE (3.5). In addition, PDE (3.6) 
simplifies to 

(1 - x 2 )(G x + GG V ) + {2p + n- 3)xG = p(p + n- 2)v + kv q . (3.11) 



We then see that the characteristic ODEs for solving this first-order PDE are 
precisely 

dv/dx — G, (1 - x 2 )dG/dx + (2p + n - 3)xG = p(p + n - 2)v + kv q , 

(3.12) 

which are satisfied due to equations (3.8) and (3.9). Hence, we have proved 
the following result. 

Lemma 1: Scaling-invariant solutions (3.7) of wave equation (2.2) correspond 
to solutions of the first-order PDE system (3.5)-(3.6) with the specific form 
(3.10). 

We now demonstrate that, in fact, a larger correspondence holds. 

Lemma 2: There is a mapping (3.1) and (3.2) between the solutions u(t,r) 
of wave equation (2.2) and the solutions G(x,v), H(x,v) of its scaling-group 
resolving system (3.5)-(3.6). The map u — > {G,H) is many-to-one and onto; 
its inverse (G, H) — > u is one-to-many and onto. 
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Proof: For any solution pair G(x, v), H(x, v), we note that the total derivative 
equation (3.3) holds by PDE (3.5), and hence G and H have the form (3.9) 
for some function u(t,r) which, however, is non-unique. Through PDE (3.6), 
we also have the total derivative equation (3.4). Substitution of expressions 
(3.9) into this equation then yields the wave equation (2.2) on u(t,r). 

To show the converse for any solution u(t,r), we express m as a function 
v(t,x) = (x/t) p u(t,t/x) through relation (3.1). Without loss of generality, we 
will assume v t is not identically zero, as otherwise we have a group-invariant 
solution v = v(x) for which Lemma 1 already establishes existence of a corre- 
sponding solution G(x, v ), H(x, v ) of the system (3.5)-(3.6). By substitution 
of the resulting expressions (3.1) into the relation (3.2), we obtain a pair 
of functions G(t,x), H(t,x) that satisfy the total derivative equations (3.3) 
and (3.4). Now we use the implicit function theorem to express t = T(v, x) for 
some function T, which exists at least locally because, by assumption, v t ^ 
in some open domain in the (t, x)-plane. Elimination of t in G and H then 
produces functions G = g(x,v) = G(T(v,x),x), H = h(x,v) = H(T(v,x),x). 
Since G and H satisfy equations (3.3) and (3.4), the pair g(x,v),h(x,v) is a 
solution of system (3.5)-(3.6). □ 

The many-to-one nature of this mapping arises when equation (3.2) is in- 
tegrated to obtain u(t,r) from G(x,v), H(x,v). Through equation (3.1), the 
integration consists of solving a consistent pair of parametric first-order ODEs 

Ut = r^Git/r, r- p u), u r = r^H^t/r, r-- p u) (3.13) 



whose general solution will involve a single arbitrary constant. 

The PDE system (3.5)-(3.6) for G(x, v ), H(x, v ) is well-suited to solve by a 
separation of variables ansatz if we exploit the feature that the nonlinearity 
in this system is only quadratic, while the interaction term given by a power 
nonlinearity u q in the wave equation now appears only as an inhomogeneous 
term v q in the system. Thus we will consider solutions of the form 

G = g(x)v a , H = h(x)v b (3.14) 



with a power dependence on v. In outline the steps to find a,b,g(x),h(x) 
proceed as follows. First the powers a, b are determined by substituting ex- 
pressions (3.14) into the PDEs (3.5) and (3.6) and balancing powers of v, in 
particular, v a , v\ v a+b ~ l in PDE (3.5), and v a , v\ v 2 ^ 1 , u 26 " 1 , v q in PDE (3.6). 
This leads to a number of cases to consider for a, b, taking into account that 
the term v q must be balanced by at least one term among the other powers. 
Equating coefficients of like powers in each case then yields an overdetermined 
system of algebraic and differential equations for g(x), h(x). Their overdeter- 
mined nature allows all these systems to be readily solved. In those cases where 
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nontrivial solutions are found, we are left with finding u(t, r) by integration of 
a parametric ODE system (3.13), which reduces to solving separable ODEs 



u 



~ a Ut = r p{1 - a) - l g(t/r), u~ h u T = r^ 1 " 6 )" 1 h(t / V) . (3.15) 



Rather than consider all possible cases exhaustively, we will now go through in 
detail only those cases for which the system of equations obtained on g(x), h(x) 
is found to be the least overdetermined, giving the most promising chance for 
nontrivial solutions to exist. As suggested by compatibility of the pair of ODEs 
(3.15), these cases are distinguished by the relation a = b. 

Case a = b = (q + l)/2: Here, PDE (3.5) involves just the power v (q+1)/2 , 
whose coefficient then yields the ODE 

xg' + h' = 0. (3.16) 



PDE (3.6) contains two powers v^ g+1 ^ 2 ,v g , and their respective coefficients 
give the ODE 

g' + (1 - n)h + xti = 0, (3.17) 



in addition to the algebraic equation 

g 2 -h 2 = 2k/(q + l). (3.18) 



Differentiation of equation (3.18) followed by substitution of equation (3.16) 
gives (xh + g)g' = and hence 

g' = or xh + g = 0. (3.19) 



For the latter possibility g = —xh, equation (3.17) reduces to nh = 0, which 
determines h — 0. But this implies g = —xh = 0, contradicting equation 
(3.18). Therefore we must have the other possibility, g' = 0. Equations (3.16) 
and (3.17) together now yield (n — l)h = and hence, for n > 1, 

h = 0. (3.20) 



Then equation (3.18) determines 

g = ±pk/{q + l). (3.21) 
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Thus, we have a solution of PDE system (3.5)-(3.6). We remark that, 
in the less interesting case n = 1, the solution generalizes to g — 
(c 2 + 2k)/ (q + 1), h — c/ \J (q + 1), with an arbitrary constant c. 

Case a = b = q: PDE (3.5) again involves a single power, yielding the ODE 

ti + xg'-g = 0. (3.22) 

Likewise, PDE (3.6) again contains two powers, leading to the ODE 

xh' -nh + g' = k (3.23) 

and the algebraic equation 

g 2 = h 2 . (3.24) 



Clearly, from equation (3.24) we have g = ±h, and then equations (3.16) 
and (3.23) reduce to an overdetermined pair of ODEs 

(x± l)ti = k + nh = h. (3.25) 



Elimination of h' gives the algebraic equation h — k/(l — n), implying h' = 0, 
which is inconsistent. Therefore, PDE system (3.5)-(3.6) has no solution in 
the present case, for any n > 1. 

Unsurprisingly, all additional cases obtained from balancing other powers can 
be shown to yield no solutions except if n = 1. 

Because separation of variables (3.14) is successful in yielding solutions of 
PDE system (3.5)-(3.6), we will consider a more general two-term ansatz 

G = gi(x)v a + g 2 (x)v, H = h 1 {x)v b + h 2 (x)v. (3.26) 



The restriction to linear powers for the additional terms here is motivated by 
the particular form of nonlinearities occurring in the PDEs. This ansatz (3.26) 
is found to yield three more solutions for n > 1. 

Proposition 5: Through ansatzes (3.14) and (3.26), the scaling-group resolving 
system (3.5)-(3.6) for n > 1 has the solutions: 



G = ±^2k/(q + l)v (q+1)/2 , H = 0, (3.27) 
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G = 0, H = (2 - n)v ± y/(2 - n)kv l ^ n ~ 2 \ q = i— ^, (3.28) 

G = ±A;(n-3)- 1 t; (n - 1)/(n - 2) , H = (2 - n)v ± G, g = (3.29) 
G = ±V^kv- 1 + x~ 1 v, H = 0, q = -3. (3.30) 

To continue, we next write down the group-resolving systems associated with 
the translation and inversion point symmetries for the wave equation (2.2), 
and present the solutions obtained via our ansatz technique. 

It should be noted that the inversion-group resolving system is equivalent to 
the translation-group resolving system upon specialization to the case of the 
conformal power q = q c , which follows from the involution property of the 
special conformal point symmetry transformation (2.29) admitted by wave 
equation (2.2). However, the complexities of the ansatz technique differs for 
these two systems, and actually as a result we will see different solutions 
emerge. 

Obvious invariants and differential invariants for the translation symmetry 
(2.3) are simply 

x — r, v = u (3.31) 



and 



G = u t , H = u r . 



(3.32) 



The resulting first-order PDE system is given by 



G x + HG V - GH V = 0, (3.33) 
GG V - HH V -H x -(n- ^x^H = kv\ (3.34) 

with independent variables x, v, and dependent variables G, H. 

Lemma 3: Solutions of the translation-group resolving system (3.33)-(3.34) of 
the form 

G = 0, H = U'(x) (3.35) 



correspond to translation-invariant solutions u = U(x) of wave equation (2.2). 
More generally, the solutions u(t,r) of the wave equation and the solutions 
G(x,v), H(x,v) of its translation-group resolving system are related by a 
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many-to-one, onto mapping u — > (G,H), and a one-to-many, onto inverse 
mapping (G,H) — > u, through equations (3.31) and (3.32). 

Explicit solutions are readily found by the same separation of variables 
ansatzes introduced in the scaling case. However, up to a change of vari- 
able, these solutions are merely the same as the ones (3.27), (3.28), (3.29) 
derived in the scaling case that have no dependence on the scaling invariant 
x = t/r. This reflects the fact that the PDE system (3.33)-(3.34) admits a 
scaling symmetry, leading to invariant solutions of the form 

G(r, u) = i*- x g(r- p u), H(r, u) = r^hir^u), p = 2/(1 - q) (3.36) 



where g, h depend only on the scaling invariant v = r p u and satisfy the 
scaling-group resolving system (3.5)-(3.6). 

For the inversion symmetry (2.5), invariants x(t, r),v(u, t, r) are given through 
the equations 

Xi nV er.2; = (t 2 + r 2 )x t + 2rtx r = 0, (3.37) 
Xi nvcr .f = (1 — n)tuv u + (t 2 + r 2 )v t + 2rtv r = 0. (3.38) 

Integration of the characteristic ODEs of these two equations yields 

x = (t 2 - r 2 )/r, v = r {n - l)/2 u. (3.39) 

Differential invariants G(u t ,u r ,u,t,r), H( t,r), with essential depen- 

dence on u t ,u r , respectively, are obtained in a similar way from solving the 
equations 

X« cr G = 0, X« cr iJ = 0, (3.40) 
where 



X^ er = X inver . + ((1 - n)u - (1 + n)tu t - 2ru r )d/du t - ((1 + n)tu r + 2ru t )d/du r 



(3.41) 



is the prolongation of the generator X inver . from (t, r, u) to (t, r, u, u t , u r ) (see 
Refs. [1,2]). This leads to the expressions 

G = r^-^dt 2 + r 2 )u t + 2rtu r + (n - l)tu), (3.42) 
H = r^/ 2 ((t 2 + r 2 )u r + 2rtu t + (n - rf—^-u). (3.43) 
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The reciprocal relations 



m = r^- n ^ 2 (t 2 - r 2 )- 2 ((t 2 + r 2 )G - 2rtH), (3.44) 

Ur = r W 2 (t 2 _ r 2 )- 2 ((t 2 + r 2 )H - 2rtG) + — — w, (3.45) 

give rise to the following first-order PDE system 

x 2 G x + GH V - HG V = 0, (3.46) 

x~ 2 (GG v - HH V ) +H X = -p{p + l)v + kv q , (3.47) 



with independent variables x, v, and dependent variables G,H, where q = 
(n + 3)/(n - 1) and p = (1 - n)/2 = 2/(1 - q). 

Lemma 4: Conformally-invariant solutions u = r^^^Uix) of wave equa- 
tion (2.2) correspond to solutions of conformal-group resolving system (3.46)- 
(3.47) of the form 

G = 0, H = -x 2 U'(x). (3.48) 

Through equations (3.39), (3.42)-(3.45), the solutions u(t, r) of the wave equa- 
tion and the solutions G(x, v), H(x, v) of its conformal-group resolving system 
are related more generally by a many-to-one, onto mapping u — > (G, H), and 
a one-to-many, onto inverse mapping (G, H) — > u. 

The previous separation of variables ansatz again readily yields explicit solu- 
tions. In contrast to the translation case, the system (3.46)-(3.47) here is able 
to yield new solutions. 

Proposition 6: Through ansatzes (3.14) and (3.26), the conformal-group re- 
solving system (3.46)-(3.47) for n > 1 has the solution 

G = ±^k(n - l)/(n + l)xv {n+l ^ {n ~ l \ H = -{n - l)xv. (3.49) 
3. 1 Group-invariant potentials and ansatzes 

We introduce an extension of our ansatz technique by taking advantage of con- 
servation laws that arise for each of the group-resolving systems. To illustrate 
the ideas, we consider the translation case first. 

It is easily seen that each of the PDEs in the translation-group resolving 
system (3.33)-(3.34) admit a conservation law 
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{H/G) v = (-1/G) X , (3.50) 
2k 

(x n -\G 2 - H 2 )) v = (2x n ~ 1 H + —x n v q ) x , (3.51) 

where, recall, x = r and v = u are the translation invariants. Associated with 
these respective conservation laws are potential variables (see Ref. [10]) given 
by 



$ x = H/G, = -1/G, (3.52) 

y x = X n -\G 2 - H 2 )/2, = x n ~ l H + -x n v q , (3.53) 

n 

with corresponding reciprocal relations 

G = -!/$„, H = -$ x /$ v , (3.54) 



and 



G = ±W 2^ 1 -«^ :E + - -^) 2 , if = x 1_n ^ - -xw 9 . (3.55) 

V n n 



The system of PDEs (3.33)-(3.34) can be rewritten in terms of each potential 
((1 - $> x 2 )/$> v ) v = -2x x - n {x n -^ x l^ v ) x + 2kv q (3.56) 



and 



(aT n tf„ - -v q )/\i2x- 1 - n ^ x + (x~ n ^ v - -vi) 2 
n V n 



= ( x~ l / \l2x- l - n m x + (a:-"^ - ^V) 2 ] . (3.57) 



X 



Note, equations (3.52) and (3.53) each provide a one-to-many, onto embedding 
of solutions G(x,v), H(x,v) of the translation-group resolving system into 
solutions $(x,v),i&(x,v) of the respective associated potential systems. 

We remark that the inhomogeneous term in the conservation law (3.51) can 
be shifted between the two sides, altering the potential \I/ by a multiple of the 
corresponding term n ^^ x n v 9+1 . Throughout, we simply fix this freedom in 
the most convenient way for solving the resulting potential systems. 
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We will now employ a separation of variables ansatz 

$ = (f)(x)v a , = ip(x)v b (3.58) 

for finding explicit solutions of each PDE (3.56) and (3.57). As before, the 
steps consist of balancing powers of v to determine a, b, followed by solving 
an overdetermined system of algebraic and differential equations for <p{x) or 
ip(x). In particular, despite the complexity of their appearance, the form of 
these PDEs becomes quite manageable to handle after the substitution (3.58) 
is carried out. This readily leads to the results 

$ = ±y/2(q + l)/k(q - I)" V 1 ^ 2 , (3.59) 
G = T^2k/(q + l)v {q+1)/2 , H = 0, (3.60) 

and 



* = -r~, x rx n v q+1 , (3.61) 

ra(n(l - q) + 1 + q) V ; 

G = ±— 1 v*y/x*(q - l) 2 + 2(n(l - q) + 1 + q)v^, 

n(l - q) + 1 + q v 

H= - q ~ 1 _ l xv q . (3.62) 

n(l-q) + l + q V 1 

We note that the second solution (3.62) is new, while the first solution (3.60) 
reproduces the solution (3.27) already obtained in Proposition 5 directly from 
the translation-group resolving system. This duplication occurs because the 
form of G, H for the $ potential using the ansatz (3.58) is precisely the same 
as the original separation of variables ansatz (3.14). 

The situation is similar in the scaling case as well as the conformal case. For the 
conformal-group resolving system (3.46)-(3.47), there are conservation laws 



(x- 2 H/G) v = (1/G) X , (3.63) 
(x- 2 (G 2 - H 2 )) v = -2{H + p(p + l)xv + kxv q ) x , (3.64) 

where x = (t 2 — r 2 )/r and v = r^ n " l ^ 2 u are the conformal invariants, with 
q = q c = (n + 3)/(n — 1) and p — (1 — n)/2. Using the potential 

ty x = (G 2 - H 2 )/x 2 , tf„ = -2H - 2x(p(p + l)v + kv q ) (3.65) 



associated with the conservation law (3.64), we find the ansatz (3.58) leads to 
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no solutions; however, via a two-term ansatz 

* = i>i(x)v b + i> 2 (x)v 2 



(3.66) 



motivated by the form of the PDE for this potential, we obtain 

m = kxv q+1 - p 2 xv 2 (3.67) 

and 

k 

V = 2——<cv q+1 - p 2 xv 2 . (3.68) 

The first result (3.67) gives a new solution 

G = ±xJp~ 2 v 2q — kv q+1 , H = —pxv + kp~ 1 xv q , q = -, p = — - — .(3.69) 

v ' n — 1 2 

The second result (3.68) reproduces the solution already obtained in Proposi- 
tion 6. This same solution is the only one found through the other potential 

<$> x = H/(x 2 G), = 1/G (3.70) 



with use of our ansatzes. 

The scaling-group resolving system (3.5)-(3.6) has conservation laws 

(G/(H + xG- pv)) v = -(1/(H + xG- pv)) x , (3.71) 
(G + xH) x - {\{H 2 - G 2 ) - pvH + -^—v q+1 ) v = (2p + n- 1)H, (3.72) 

with x = t/r and v = r~ p u, where p — 2/(1 — q). For strict conservation 
in equation (3.72), we need 2p + n — 1 =0, which is seen to imply that q = 
(n + 3)/(n — 1) is restricted to the conformal power q c . (Note, this parallels the 
condition for scaling invariance of the action principle for the wave equation.) 
From the potential 

V v = G + xH, y x = -(H 2 -G 2 )-pvH + -^—v q+1 (3.73) 

2 q ~\~ 1 

associated with the resulting conservation law (3.72), we find by ansatz (3.58) 

= kq^l q+1 ( 4) 

2g+l V ; 
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and this gives a solution 



G = p 




H = --(q-l)vi + x- 1 G, 



Q = 



n + 3 



1 — n 
2 



(3.75) 
(3.76) 



Our results are summarized as follows. 

Proposition 7: By means of potentials introduced through conservation laws, 
the group-resolving systems associated with translations, scalings, and inver- 
sions have the additional solutions (3.62), (3.76), (3.69), for n > 1. 

3.2 Optimal group-resolving systems 

So far we have considered the group foliation resolving systems that arise from 
the geometrical subgroups of scaling symmetries, translation symmetries, and 
(in the conformally-invariant case q = q c ) inversion symmetries admitted by 
wave equation (2.2). In the case of the conformal power, there is an essen- 
tially different group foliation and associated resolving system given by a sub- 
group consisting of translation symmetries combined with inversion symme- 
tries. Note, in general the criterion for foliations to be equivalent with respect 
to the full admitted symmetry group is that they must come from subgroups 
related by conjugation. For the group-theoretical reasons discussed earlier in 
connection with classical invariant solutions (cf. Sec. 2.4), an enumeration of 
conjugacy classes under the full group of geometrical point symmetries of the 
wave equation (2.2) is represented by the symmetry generators X sca i , X trans ., 
in addition to X invcr + X trans . in the conformally-invariant case. 

To proceed in the latter case, the invariants and differential invariants of the 
point symmetry X invcr . + X trans . are straightforward generalizations of those 
derived for inversions. We have 



x = 



(l + t 2 -r 2 )/r, v = r- p u, 



(3.77) 



and 



G = r~ p ((l + t 2 + r 2 )u t + 2rtu r - 2ptu), 



(3.78) 
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H = r- p ((l + t 2 + r 2 )u r + 2rtu t - pr _1 (l + t 2 + r 2 )w), (3.79) 

where p — (1 — n)/2 — 2/(1 — g) and q = (n + 3)/(n — 1) = g c . Through the 
reciprocal relations 

^ = r p (l + (t 2 - r 2 ) 2 + 2(t 2 + r 2 ))- 1 ((l + t 2 + r 2 )G - 2r*J^80) 
Ur = r p {\ + (t 2 - r 2 ) 2 + 2(t 2 + r 2 ))- 1 ((l + t 2 + r 2 )i7 - 2rtG) + pr^l) 

we obtain the first-order PDE system 

(4 + x) 2 G x + GH V - HG V = 0, (3.82) 
(4 + x 2 )-\GG v - HH V ) + H X = -p(p + l)v + kv q , (3.83) 

which differs from the conformal-group resolving system (3.46)-(3.47) only by 
the factor x 2 + 4 in place of x 2 . 

As before, we first try separation of variables ansatzes (3.14) and (3.26) to solve 
the group-resolving system (3.82)-(3.83). However, unlike the conformal group 
case, this yields no solutions here. Next we try exploiting the conservation law 
admitted by each of the PDEs (3.82) and (3.83), 

{{x 2 + A)- 1 H/G) v = (1/G) X , (3.84) 
((:r 2 + 4)- 1 (G 2 - H 2 )) v = -2(H + p(p + l)xv + kxv q ) x . (3.85) 

Through the potential $ x = H/(x 2 G), & v = 1/G, associated with con- 
servation law (3.84), no solutions are found via the ansatz (3.58) or its 
two-term generalization. Using the potential ty x = (G 2 — H 2 )/(x 2 + 4), 
ty v = —2H — 2x(p(p + l)v + kv g ), from conservation law (3.85), finally we 
find a two-term ansatz (3.66) leads to solutions analogously to the case of the 
conformal-group resolving potential system. 

Proposition 8: Group-resolving system (3.82)-(3.83) has the solutions for n > 
1 

$ = kxv q+l - p 2 xv%3.86) 
G = ±x\J (4 + kp~ 2 x 2 v q ~ l )[kxv q+1 — p 2 x 2 v 2 ), H = pxv — kp~ 1 xv\$.87) 

and 

2k 

* = -—^xv q+1 - p 2 xv 2 , (3.88) 
G = ±x^2k(q + l)-!(4 + x 2 )vi +l - 4p 2 v 2 , H = pxv. (3.89) 
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4 Exact solutions 



To construct solutions u(t,r) of wave equation (2.2) from solutions G(x,v), 
H(x,v) of its symmetry group-resolving systems, we integrate the corre- 
sponding group-invariant equations that define the symmetry invariants x, v, 
and differential invariants G, H used in formulating those systems in terms 
of t,r,u,u t ,u r . These group-invariant equations (3.2), (3.32), (3.44)-(3.45), 
(3.80)-(3.81) take the form of linear, inhomogeneous systems of two first- 
order PDEs, in which G(x,v), H(x,v) are explicit expressions. Because of the 
inherent symmetry possessed by such systems, they reduce to a parametric 
pair of first-order ODEs whose integration to obtain u(t, r) is straightforward. 
We will now list all the explicit solution expressions found for u(t,r). 

Theorem 1: Wave equation (2.2) has the following exact solutions arising from 
the explicit solutions of its group-resolving systems obtained in Propositions 5 
to 8 for n > 1: 



2/(1-9) 



-=(±^^(.-l)(t + c)j , ^-1, (4.1) 

u=\ , , . ' -((t + c) 2 -r 2 )) , q^^—, n+ 1,(4.2) 

\2{q{\ -n) +n + 1) VV ; ') ' y n-l' T ' V ; 

n-3 \^ 2 ^ 4 _ n 



u=[±V-k^-^ s r + cr 3 ~ n ) , q=^, n^2,3 (4.3) 

\ \/n — 2 / n — z 



u = V±2V-A;*(l + ct), u = -yj±2y/-kt(l + ct), q = -3, (4.5) 

(l-n)/2 



u = 



(/I / 1 \ n )/^ Q 



where c is an arbitrary constant. 

Solutions (4.5) to (4.8) can be slightly generalized by applying a time- 
translation to t. It is easy to verify that all the resulting families of solutions, 
with one exception, are invariant with respect to the full geometrical symme- 
try group admitted by the wave equation (2.2), depending on the nonlinearity 
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power q. This group comprises translations, scalings, and inversions for solu- 
tions (4.6) and (4.7) involving the conformal power q c = (n + 3)/(n — 1), and 
otherwise only translations and scalings for solutions (4.1) to (4.5) involving 
other powers q ^ q c . The solution (4.8) is exceptional as it fails to be invariant 
under inversions, despite falling within the conformal power case. We note 
these invariance properties are determined by the group-invariant variables of 
the PDE systems from which the individual solutions arise. In particular, the 
lack of conformal invariance for solution (4.8) is accounted for because it comes 
from the scaling-group resolving system rather than the conformal-group re- 
solving system. Thus, this solution can be generalized to a two-parameter 
family under inversion transformations (A. 3), yielding 



Ah 1 9 \(i-«)/4 

u = f ( ± r 2 + \ (c(t 2 - r 2 ) ±c-\l + Tct + 5 2 (t 2 - r 2 ))) 2 ) ^ 



with arbitrary constants c, c. 

We remark that, once translations and scalings are taken into account, the 
solution families (4.6) and (4.7) coming from the conformal-group resolving 
system in Propositions 6 and 7 are the same as the ones obtained through the 
translation-inversion group resolving system in Proposition 8. 

For comparison with classical symmetry reduction, note that solution (4.3) is 
of translation-invariant form, and solutions (4.1), (4.2), (4.4) are of scaling- 
invariant form to within a time translation, as also are solutions (4.6) and (4.7) 
to within an inversion transformation. Of interest is that the two solutions (4.5) 
and (4.8) do not have such group-invariant forms with respect to any of the 
geometrical point symmetries of wave equation (2.2) (namely, under the action 
of the full group of geometrical symmetries, they do not lie on the orbits of 
any scaling-invariant, translation-invariant, or inversion-invariant solutions). 



4-1 Analytical behavior 



It is certainly noteworthy here that we have found explicit solutions (4.3) to 
(4.5) for three special nonlinearity powers, q = q = q = —3, of 
which the last two are not distinguished by any special cases of the symmetry 
structure of the wave equation or of the ODEs for group-invariant solutions, 
for n > 1. All these powers are subcritical, i.e. q < g* where g* = (n + 
2)/(n — 2) is the critical power. Two other solutions (4.1) and (4.2) hold for 
almost all powers q, including in particular the critical power as well as all 
supercritical powers. However, the analytical behavior of these two solutions 
is not so noteworthy, as u blows up at t = in solution (4.1) and at the light 
cone \t\ = r in solution (4.2). Related to this blow up, the initial data of both 



25 



solutions at earlier times, t = t < 0, has infinite energy whenever q > g*, as 
seen from the energy integral 



°° 1 1 

£=f (^ 2 + ~ 



q+l 



k 



u 




dr. 



(4.10) 







Note S will be conserved, dS/dt = 0, for any solutions u(t,r) with smooth 
initial data having suitable asymptotic decay for large r. 

The analytical features of the solutions with subcritical powers are worth com- 
menting on. First, solution (4.5) is radially homogeneous (i.e. r independent) 
and hence it satisfies the nonlinear oscillator ODE uu = ku~ 3 . There is no 
blow up of u, but u t becomes singular at t = as well as t = — c _1 in this 
solution, which is caused by the negative nonlinearity power. Next, solution 
(4.3) is static (i.e. t independent). Asymptotically it is badly behaved for large 
r since u does not decay. Also it is singular at r = for c 7^ 0, so this solution 
thus describes a monopole with infinite energy. Note these two solutions are 
real- valued in the case k = — 1, precisely when the energy integral (4.10) is 
positive definite. The remaining non-static solutions (4.4) to (4.8) with sub- 
critical powers q < q* all exhibit a more interesting blow up of u. 

We recall that the blow-up set of a solution u(t,r) is defined to be those 
points in the (r, t) half-plane where 11 or a derivative of u is singular. For 
solutions (4.4) and (4.7), the blow-up sets are given by the light cone r = \t\, 
in addition to, respectively, the timelike line r = and the shifted light cone 
r — \t + c _1 |, c ^ 0. The blow-up set for solution (4.8) in the "+" case is 
a pair of shifted light cones r = \t ± c _1 |, c ^ 0. Because all these blow-up 
sets contain a light cone, the solutions describe initial data with a singularity 
at r = I to + c_1 | > that propagates forward with unit speed in the (r,t) 
half-plane. This propagation reflects the fact that the wave equation (2.2) is a 
semilinear hyperbolic PDE with its principal (second derivative) part given by 
the radial wave operator d\ — <9 2 whose characteristic curves are light cones. 

Of the most interest analytically is solution (4.6) with the conformal power 
q = q c . For k — +1 and c < this solution is real- valued. Its blow-up set is a 
hyperbola \t + 1*\ = \/r 2 + 1* 2 , where t* = ^^-|c| _1 , which is seen to have a 
timelike future branch starting at t — and a spacelike past branch ending at 
t = —2t*. At intermediate times > t > —2t*, there is no blow up for all r > 0, 
and so u is globally smooth in this strip in the (r, t) half-plane. Consequently, 
for any fixed time t = t , with > t > —2t*, we have smooth initial data 
that evolves to a blow-up at time t — 0, as illustrated in Figure 1. This initial 
data, moreover, is found to have finite energy, due to the asymptotic decay 
u = | c |( 1 - n )/ 2 r 1 - n -f 0(r _n ) for large r. But, because of the sign of k = +1, the 
energy associated with solutions is not positive definite, which is what allows 
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Blow-up set and light cones 




Fig. 1. 

the blow-up to occur. 

Thus, this solution (4.6) describes a finite-time blow up for the conformal wave 
equation 

u tt - u„ - (n - l)u r /r = fc u 1+2/(n-1) , k = +1 (4.11) 
whose nonlinearity is unstable due to its sign. The blow-up rate is 

u = (2| c |O (1 ~ f0/2 * (1 ~ n)/2 + 0(r (1+n)/2 ) 

for t — > at r = for solution (4.6). A worthwhile question is to what extent 
this behavior may be an attractor to which other blow-up solutions evolve as 
the blow up time is approached. 

Also worthy of analytical interest, the solution (4.8) in the "— " case has no 
blow up anywhere in the (r, t) half-plane. For large r it behaves like 

u = (Vk\c\/(n - i))(i-")/V-" + 0(r- n ) 

and this asymptotic decay is sufficient to ensure that its energy is finite. More- 
over, the initial data for this solution is smooth for k = +1 and evolves with 
a long-time behavior given by 

u = (2Vkc/(n - + 0(r n ) 

for large t. It is thus a global smooth solution of finite energy for the conformal 
wave equation (4.11). 
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The contrasting behaviors of the previous two solutions with k = +1 for 
the conformal power q = q c is typical of wave equations with an unstable 
nonlinearity [7]. 

We remark that solution (4.9) has a similar analytical behavior to solution 
(4.8). 



4-2 Group-invariant solutions 



Finally, we will also consider for comparison the solutions (2.13), (2.16), (2.21) 
derived by classical symmetry reduction of the wave equation (2.2). These so- 
lutions were found to satisfy nonlinear, unstable oscillator equations, involving 
the special powers q = q = f° r the nonlinearity. It is possible to carry 
out the quadratures to obtain the solutions in explicit form in the special case 
when the integration constant corresponding to the oscillator energy is zero. In 
terms of the canonical coordinates x, v adapted to the symmetry reductions, 
the quadratures (2.21) and (2.13) both yield 



,= l^^(l-(tanhx) ±2 l , (4.12) 



(n-2)/4 



while the other quadrature (2.16) gives 



. exp(± v /z A;x) , for n = 3, 

[±yj^(n-3)x) , for 2, 3. 



Hence we obtain scaling-invariant solutions 



and translation-invariant solutions 



n=(±^iy~ 2) '\r^ir^, 8 =2±|, n + % (4.16) 

(n-2)/(n-3) 
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(we omit the case n = 3 for the latter solution as the power q — 1 is trivial). In 
addition, we also list the special invariant solutions (2.14) and (2.26) arising 
from the extra point symmetries admitted by the translation, scaling, and 
inversion symmetry reductions. These yield 



(2^) 



(n-2)/2 



n-2 



u 



^-l)(^-3) V W " 1)/4 ^ /2 



4k 



n + 3 
n — 1 



(4.18) 
(4.19) 



We note solutions (4.14), (4.15), (4.17) are specializations of solutions (4.1) 
to (4.3). Thus, the only extra solutions arising from symmetry reduction are 
the ones (4.16), (4.18) and (4.19), all of which are static. 

The two static solutions (4.18) and (4.19) clearly are singular at r = and 
hence describe monopoles of wave equation (2.2) with the critical power q = q* 
and conformal power q = q c , respectively. However, it is readily seen that their 
energy is infinite due to the asymptotic rate of decay of u with r. 

In contrast, the other static solution (4.16) with critical power q = q* has 
finite energy and is non-singular at r = 0. In the case k = — 1, where the 
nonlinearity of the wave equation is stable, this solution is singular at radius 
r = 1 (note this location can be moved to any non-zero radius r = r* by 
applying a scaling symmetry to the solution). In the opposite case k — +1, 
there is no singularity at any radius r > and the solution is globally smooth, 
but here the nonlinearity is unstable. Thus, this solution (4.16) describes a 
finite-energy static "soliton" for the unstable critical-power wave equation 

u tt -u rr -{n-l)u r /r = ku 1+Al{n - 2 \ k = +1. (4.20) 



The existence of such a solution suggests that there may be very interesting 
analytical long-time behavior for smooth finite-energy solutions of this wave 
equation. 



5 Remarks 



Group foliations in general convert a nonlinear PDE into an equivalent sys- 
tem of first-order equations (the group-resolving system) by means of group- 
invariant variables associated with an admitted group of transformations. The 
construction of group foliations using admitted point symmetry groups of non- 
linear PDEs is originally due to Lie and Vessiot and was revived in modern 
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form by Ovsiannikov [4]. Because the equations for a group foliation contain 
all solutions of the given nonlinear PDE, ansatzes or differential constraints 
must be used to reduce the equations into an overdetermined group-resolving 
system to obtain explicit solutions. Compared with classical symmetry re- 
duction, a main difficulty to-date has been how to find effective, systematic 
ansatzes that lead to useful reductions. 

Recent work [6,11] has introduced an ansatz based on imposing commutativity 
of operators for group-invariant differentiation, developed for the example of 
the infinite-dimensional conformal symmetry group of the "heavenly equation" 
that describes the geometry of self-dual gravitational fields in General Rela- 
tivity. This ansatz is shown to give a nontrivial reduction of the conformal- 
group foliation equations, yielding explicit exact solutions. Other work, on 
hydrodynamic equations with infinite-dimensional symmetry groups, appears 
in Ref. [4,12]. 

The semilinear n-dimensional radial wave equation we have considered is an 
example with a finite-dimensional symmetry group. Our results show that a 
basic separation of variables ansatz is effective in reducing the group foliation 
equations, especially when combined with the use of potentials based on con- 
servation laws of these equations. It is significant that we find some reductions 
give exact solutions having interesting analytical properties. 

For future work, applications of these methods to wave equations with non- 
linear wave speeds and to nonlinear heat conduction equations, as well as to 
integrable nonlinear evolutionary equations would be fruitful to consider. It 
is expected that in general our method should be well suited to equations 
with power nonlinearities. Concerning group-theoretical aspects, it would also 
be worthwhile to examine ansatzes for reduction of group-resolving systems 
in connection with the well-known non-classical method [13,14,15] for finding 
exact solutions of nonlinear PDEs through invariant surface conditions. 
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A Appendix 



Here we state the form of the translation subgroup, scaling subgroup, and in- 
version subgroup comprising the one-dimensional geometric point symmetry 
groups admitted by wave equation (2.2) for nonlinearity powers q ^ 1. Their 
action on solutions is described by the following one-parameter (A) transfor- 
mations: 



translations u(t, r) — > u(t + A, r) for all q, (A.l) 
scalings u(t,r) — > \~ p u(\t, Ar) for all q, (A. 2) 

inversions u(t, r) — > f^^w^A)^ + \(t 2 — r 2 ), f2(A)?") only for q (A<?3,) 



with Q 



t 2 - r 2 



(A) 



(t + X(t 2 - r 2 )) 2 - r 



2 ' 



where p = 2/(1 — q). We note that, in terms of the invariants £ = t/r, x = 
(t 2 -r 2 )/r, v = r p u associated with scalings and inversions, the one-parameter 
inversion transformations (A. 3) act in a simple way by 

v(£,x) ^v(£ + \x,x) (A.4) 



where v(£,x) = r~ p u(t,r) for solutions u(t,r). Finally, if we combine an in- 
version (A. 3) for A = 1/e and a scaling (A. 2) for A = e 2 together with a 
translation (A.l) for A = e, the limit e — > yields the special conformal trans- 
formation (2.29) belonging to the full point symmetry group SX(2,R) in the 
conformal power case q = q c . This conformal symmetry is distinguished by 
providing an involution under which the translation and inversion subgroups 
are exchanged while the scaling subgroup is kept invariant. 
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